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Abstract 

Classical electrodynamics can be based on the conservation laws 
of electric charge and magnetic flux. Both laws are independent of the 
metric and the linear connection of spacctimc. Within the framework 
of such a premetric electrodynamics — provided a local and linear 
constitutive law of the vacuum is added — the propagation of elec- 
tromagnetic waves in the geometric-optics limit can be studied. The 
wave vectors of the wave fronts obey a quartic extended Fresnel equa- 
tion. If one forbids birefringence in vacuum, the light cone emerges 
and Maxwell-Lorentz vacuum electrodynamics can be recovered. If 
minimal coupling of electrodynamics to gravity is assumed, then only 
the gravitational potential, i.e., the metric of spacetime, emerges in 
the constitutive law. We discuss recent results within this general 
framework. 
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1 Introduction 

Premetric classical electrodynamics has been consistently formulated in 
2003 in the book of Hehl and Obukhov 2 . Related developments can 
be found in the recent books of Lindell [3] and Russer [3] and the articles 
of Delphenich [31 [f)| • One of the basic ideas is to formulate electrodynamics 
in such a way that the metric of spacetime (that is, the gravitational po- 
tential of Einstein's theory of gravity) doesn't enter the fundamental laws 
of electrodynamics. The practicability of such an approach has been shown 
in detail in [2J , see also the earlier work of Toupin & Truesdell [S] , Post 
Kiehn et al. and Kovetz |12j . In this lecture, we concentrate 
on the post 2003 development, for earlier results we refer to the extended 
bibliography in [2]. 

After the publication of [2J , premetric electrodynamics turned out to be 
a lively and active subject. Mainly the following aspects were developed 
(we order them roughly chronologically): 

1. Inclusion of magnetic charges in a consistent and metric- free way by 
Kaiser |13j and by Hehl and Obukhov |14| . earlier work had been 
done by Edelen [To] . 

2. Application of the premetric formalism to the Quantum Hall Effect 
and prediction of the independence of the QHE on an external grav- 
itational by Hehl, Obukhov, and Rosenow |16j . 

3. Possible non-minimal coupling of gravity to electromagnetism , in par- 
ticular to the torsion field, by Solanki, Preuss, and Haugan [171 ITS] , 
see also Hehl and Obukhov j!9j, an d, subsequently, b y using the for- 
malism of (JJ, by Rubilar et al. [20] and by Itin et al. |21L 122 j . 

4. Derivation of the signature of the metric of spacetime from the Lenz 
rule (connected with Faraday's induction law) by Itin et al. [23] . 

5. The Lo rent z- violating CF J -electrodynamics, see Carroll, Field, and 
Jackiw [24j , was reformulated in the premetric formalism by Itin [25j . 
It enhances the transparency of the CFJ-approach. 

6. New derivation by Lammerzahl et al. J2JJ, see also |27J, of the metric 
from linear premetric electrodynamics: We forbid birefringence in 
vacuum. 

7. New insight in the physical dimensions of electrodynamics and in 
the time variability of the vacuum impedance and the fine structure 
constant by Tobar [2S] and Hehl & Obukhov 
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8. Effect of the skewon field on light propagation in an quite detailed 
article by Obukhov et al. |5D| . 

9. Relation between axion electrodynamics, see Wil czek |31| . and the 
so-called Post constraint, see Hehl and Obukhov [32]; compare the 
perfect electromagnetic conductor (PEMC) of Lindell and Sihvola 
|33| I34j . Scattering of electromagnetic fields at the PEMC. 

10. Symbol and hyperbolic polynomials of the wave equation and its re- 
lation to premetric electrodynamics as discussed by Beig [3*5] . 

1 1 . Nonlinear electrodynamics by Delphenich > see a ls° |37l I38j . 

We will now integrate these topics into a short systematic review of 
premetric electrodynamics. 

2 Axiom 1: Electric charge conservation 

We base electrodynamics on electric charge conservation. The conserva- 
tion of electric charge was already recognized as fundamental law during 
the time of Franklin (around 1750) well before Coulomb discovered his force 
law in 1785. Nowadays, at a time, at which one can catch single electrons 
and single protons in traps and can count them individually, we are more 
sure than ever that electric charge conservation is a valid fundamental law 
of nature. Therefore matter carries as a primary quality something called 
electric charge which only occurs in positive or negative units of an elemen- 
tary charge e (or, in the case of quarks, in l/3th of it) and which can be 
counted in principle. Thus it is justified to introduce the physical dimension 
of charge q as a new and independent concept. 

Mathematically, the first axiom is most conveniently formulated in 
terms of exterior calculus with the help of the 3-form J of electrical cur- 
rent. Its integral over an arbitrary 3-dimensional domain describes the total 
charge contained therein. It can be determined by counting the particles 
carrying an elementary charge. The latter can be taken as a unit of charge, 
in the sense of the theory of dimensions. The conservation of the current 
reads, in exterior calculus, 

dJ = 0. (1) 

This law is metric-independent since it is based on a counting procedure for 
elementary charges. Given a spacetime foliation, we can 1+3 decompose the 
current J = —j Ada + p into the 2-form of the electric current density j and 
the 3-form p of the electric charge density. Then the charge conservation 
law assumes a more familiar form of the continuity equation 

p + dj = 0, (2) 

where the dot denotes the time derivative (Lie derivative with respect to 
a suitable vector field) and d the 3-dimensional exterior derivative. Both 
(^Q) and © can be equivalently formulated in the integral form, see j2] for 
details. 
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Because of the first axiom and according to a theorem of de Rham, 
we can represent the conserved electric current as an exterior derivative (a 
generalized "divergence"), i.e., 

J = dH. (3) 

In this way, we naturally introduce the 2-form H of the electromagnetic 
excitation. Certainly, this quantity is not uniquely determined from 
since the transformation H — > H + dtp evidently preserves the form of 
However, this arbitrariness is eventually fixed in the actual measurements 
performed with the help of ideal conductors and superconductors, which 
means that the electromagnetic excitation does have a direct operational 
significance, see Raith |3*5] . 

The (1 + 3)-decomposition of H is obtained similarly to the decompo- 
sition of the current J, thereby introducing the 2-form T> of the electric 
excitation (historical name: "dielectric displacement") and the 1-form 7i 
of the magnetic excitation ("magnetic field"): 

H = —H Ada + V. (4) 

Inserting © as well as the 1 + 3 decomposition of the current into (j3J), we 
recover the pair of the 3-dimensional inhomogeneous Maxwell equations 

dH = J { . iV=P > (5) 
\-V + dH =j. 

Charge conservation is a law that is valid in macro- as well as in micro- 
physics. If charge conservation was violated, we could introduce a 4-form 
N (one component only!) according to dJ = N. Lammerzahl, Macias, 
and Miiller @U] studied a seemingly more restricted class of models vio- 
lating charge conservation. Such models can be used as test theories for 
experiments on the checking of charge conservation. 



3 Axiom 2: Lorentz force density 

With charge conservation alone, we arrived at the inhomogeneous Maxwell 
equations © • Now we need some more input for deriving the homogeneous 
Maxwell equations. In the purely electric case with a test charge q, we have 
in terms of components 

F a ~qE a , (6) 

with F as force and E as electric field covector. Generalizing (jHJ), the 
simplest relativistic ansatz for defining the electromagnetic field reads: 

force density ~ field strength x charge current density . (7) 

We know from conventional Lagrangian mechanics that the force ~ dL / dx l 
is represented by a covector with the dimension of h/[dx l ] = (h/t, h/i) (here 
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h, t, and I denote the dimensions of action, time, and length, respectively). 
@@@ We know from Lagrangian mechanics that the force ~ dL/dx 1 is 
represented by a covector with the absolute dimension of action h (here 
h is not the Planck constant but rather only denotes its dimension) M@@ 
Accordingly, with the covectorial force density f a , the ansatz (J7J) can be 
made more precise. The Lorentz force density acting on a charge density 
is encoded into Axiom 2: 

f a = (e a \F)AJ. (8) 

Here e a is an arbitrary, i.e., anholonomic frame or tetrad, a basis of the 
tangent space, with a = 0, 1, 2, 3, and J denotes the interior product (con- 
traction). The axiom (jHJ) should be read as an operational procedure for 
defining the electromagnetic field strength 2-form F in terms of the force 
density f a , known from mechanics, and the current density J, known from 
charge conservation. The 1 + 3 decomposition of F reads 

F = EAda + B, (9) 

thus defining the electric and the magnetic field strengths E and B (a 
differential 1-form and the 2-form, respectively). 



4 Axiom 3: Magnetic flux conservation 

Let us recall that the (four dimensional) electric current 3-form J decom- 
poses as J = —j A da + p. Accordingly, f c p represents the electric charge, 

p the electric charge density, and j the corresponding current density in 
three dimensions. Analogously — apart from the fact that J is a 3-form 
and F a 2-form — J^, B represents the magnetic flux and B the magnetic 

flux density. Consequently, E has to be interpreted as the current of the 
magnetic flux density. We could call F the (four dimensional) current of 
magnetic flux. 

Why isn't that conventionally presented in textbooks in this way? There 
seem to be two reasons: (i) In textbooks the 4-dimensional point of view is 
only a supplementary structure which is used eventually in order to crown 
the development of electromagnetism, whereas we take, in the approach 
presented, the 4-dimensional view on electromagnetism as fundamental and 
indispensable right from the beginning, (ii) A "charge" density is assumed 
intuitively to be a 3-form, i.e., one has to integrate the charge density over 
a 3- volume in order to find a net charge. This is a widespread prejudice. 
However, also a 2-form like B, if integrated over a 2-dimensional surface, can 
be considered to be a "charge" (here the magnetic flux) and, accordingly, 
a 1-form like E as the corresponding "current density" (here the current of 
the magnetic flux density). Therefore, Axiom 1 in local form, dJ = 0, as 
electric charge conservation, has as analog — in local form — 



dF = { . ~ ' (10) 
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the law of magnetic flux conservation. Since t he Faraday induction law 
and the sourcelessness of B are consequences of (|10|) . this axiom has a firm 
experimental underpinning. Note that the induction law has the form of a 
continuity equation for the "charge density" B. 

Nevertheless, similar as with the possible violation of electric charge 
conservation, there have been numerous attempts to postulate the existence 
of magnetic charge according to 



with the magnetic current density K = — J m g A da + p mg . 

It has bee n sh own by Edelen IT51 . and more completely by Kaiser |13j 
and two of us that the law can be made to fit into the premetric 
approach, provided electric charge conservation is fulfilled. 

Clearly, if besides the electric charge p c \, a magnetic charge p mg is 
present, we have to expect that the electromagnetic field exerts a force 
on it. The absolute dimension [J] x [F], which enters the Lorentz force 
density (JHJ), is electric charge x magnetic flux = action. The same dimen- 
sion can be found for the expression [K] x [H] , namely [K] = magnetic flux 
and [H] = electric charge. We recognize here a certain reciprocity between 
electricity and magnetisms. The revised Axiom 2, taking into account the 
linearity of electrodynamics, would then read 



The necessity of the minus sign will be explained in the next section. There 
it will also be shown that the expression for the energy-momentum current 
of the electromagnetic field remains the same, with or without magnetic 
charge, a quite remarkable fact that enables us to revise the premetric ax- 
iomatics in such a way as to accommodate the possible existence of magnetic 
charge. 

However, we shouldn't forget that all experiments done so far attest to 
the absence of magnetic charge. Hence the experimental situation leaves 
no doubt about the absence of magnetic charge in nature. Even more so, 
theoretical re ason s make such a structure also implausible: The modified 
Axiom 2, see (|12|). could no longer be read as a clear cut operational def- 
inition of the electromagnetic field strength F. It would mix with H in a 
rather confusing way. 

In the 19th and the first half of the 20th century, in spite of Ampere's 
hypothesis of the electric origin of magnetic effects — which today is no 
longer a hypothesis but rather an experimental fact — magnetostatics was 
built up in analogy to electrostatics. The force on the fictitious magnetic 
charges were then given by p mg A TL. This coincides with the leading piece 
of the last term of the right-hand-side of 1)12(1 . This is the reason why 
TC was called the magnetic field strength in analogy to the electric field 
strength E, with p e \ AE as Lorentz force density, even though the magnetic 
charges were assumed to be fictitious. We recognize here once more that 
the premetric formalism induces transparency into electrodynamics. 



dF = K 




(11) 




(12) 
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5 Axiom 4: Energy-momentum current 

The field equations of electrodynamics dH = J and dF = have been 
set up. Now we have to turn to the energy- momentum distributio n in 
the electromagnetic field. Since the Lorentz force density (JSJ) [or (|12|)] 
determines the relation between the electromagnetic field and mechanics, 
this formula must be the starting point for a discuss ion of energy and 
momentum. If we introduce for the magnetic piece in 1)12(1 a factor /3, to 
be determined later, then 

f a = (e a \F) A J-f3(e a \H) AK. (13) 

For K = 0, we recover the case of the absence of magnetic charges . W e 
assume magnetic charges and substitute the Maxwell equations into (|13|) : 

f a = (e a \F)AdH-(3(e a \H)AdF. (14) 

Note, however, that this equation is also valid for ordinary electrodynamics 
(withou t m agnetic charges) since under those conditions dF = 0. In other 
words, (|14|) is valid in both cases, with or without magnetic charges. The 
same will be true for the subsequent for mula s. 

We integrate partially both terms in (|14|) . see 0: 

f a = d[f3FA(e a \H)-HA{e a \F)] - (3 F A d{e a \ H) + H A d(e a \ F) . (15) 

The expression under the exterior derivative has already the desired form. 
We recall the main formula for the Lie derivative of a n arbitrary form <£, 
namely £ ea <S> = d(e a J$) + e a \(d$) , see Frankel |4"T| This allows us to 
transform the second part on the right-hand-side of (|15|) : 

f a = d[(3F A{e a \H) - HA(e a \F)] 
-(3FA(£ ea H) + HA(£ ea F) 
+(3FAe a \{dH) - H A e a \ (dF) . (16) 

The last line can be put into the form 

+/3 e a \ [F AdH]-f3 (e a \ F) A dH — e a \ [H A dF] + (e a \ H) A dF . (17) 

The expressions in the square brackets are 5-forms and vanish. Two terms 
are left over, and we find 

f a = d[(3F A{e a \H) - HA(e a \F)] 
-(3FA{£ ea H) + HA{£ ea F) 
-P{e a \F)AdH + (e a \H)AdF. (18) 

With the help of (|14|) . the last line can be rewritten as 

-pf a + (l-p 2 )(e a \H) AdF. (19) 
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Thus, (TEJ) reads 

{l + (3)f a = d[(3F A(e a \H) - HA(e a \F)] 

—(3 F A (£e a H) + H A (£ eoi F) 

+{1- (3 2 ){e a \H) A dF. (20) 

As we mentioned above, also this formula is correct with and without 
magnetic charges as long as electric charge conservation dJ = is taken 
for granted. Without magnetic charge the last term drops out because 
of dF = 0. Then the choice of (3 = 1 yields the conventional energy- 
momentum current of Minkowski. With magnetic charge the last term 
disturbs the whole set up. We cannot find a reasonable energy- momentum 
curre nt in this case unless we put f3 = ±1. However, f3 = — 1 would trivi- 
alize (|18|) to a mathematical identity and the Lorentz force density would 
be lost. Hence ( 5 = +1 is the only reasonable option. In other wo rds, we 
choose (3 = 1 in (|13|) thus arriving at what we displayed already in (|12|) . 
Consequently, in both cases, we arrive at the same relation 

f a = d k Z a + X a , (21) 
with the kinematic energy-momentum 3-form of the electromagnetic field, 

k S a := i [F A (e a \ H) - H A (e a \ F)\ (22) 

and the force density 4-form 

X a :=-~(FA£ ea H-HA£ ea F). (23) 

For the deri vation of the energy-momentum current, we could a ltern atively 
require (|21|) right from the beginning. Then the first line in 1)20(1 is rea- 
sonable because of the exterior derivative, the second line has to emerge 
because of the Lie derivative (in a special frame it is zero). Accordingly, 
the third line has to be zero. 

The energy- momentum current (|22|) in premetric exterior calculus ha s 
also been derived by Lindell & Jancewicz [32]) Segev [33], and Kaiser [T5] . 
see also the modified derivation by Itin et al. |21[ I23j . 

The energy- momentum localization, as specified in (|22j) . represents our 

Axiom 4. The current k S Q is electric/magnetic reciprocal under the sub- 
stitutions H — > (F , F — > —^H, with an arbitrary pseudo-scalar function 
as discussed in detail in [2]. There it is also shown that X a vanishes in 
specific cases making then k S a to a conserved quantity. 

6 Premetric electrodynamics and the Maxwell 
equations 

Until now, all our considerations are generally covariant and metric-free 
and connection-free. They are valid in flat Minkowskian and in curved 
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Riemannian spacetime, that is, in special relativity theory (SR) and in 
general relativity theory (GR), and even in a spacetime possibly carrying 
torsion and/ or nonmetricity. Therefore Maxwell's equations in the form of 
((SJ) and (|1U|) represent the optimal formulation of the fundamental laws of 
classical electrodynamics. Nonminimal coupling s, w hich induce additional 
fundamental constants, will be discussed in Sec ll2l In order to complete 
the theory, we need a relation between H and F, which we will discuss in 
the next section. However, before that we will have a look at dimensional 
analysis, today a largely neglected subject. 

In accord ance with the general dimensional analysis of Schouten and 
Dorgelo, see |44j . and, in particular, of Post [S], the absolute dimension 
of J is that of a charge: [J] = q (since the integral of the current over a 
3-dimensional spatial domain yields the total charge). Since the exterior 
derivative is dimensionless, [d] = 1, and since the electromagnetic excitation 
is given by the inhomogeneous Maxwell equation ©, we conclude that the 
absolute dimension of the excitation [H] = q. Then, by means of @, the 
absolute dimensions of the electric and magnetic excitations turn out to 
be [T>] = q and [TC] = q/t. The relative dimensions are those of their 
frame components, [7i a ] = q/(t£) and [V^] = q/l 2 , with the spatial indices 
a,b = 1,2,3. These are the "physical" dimensions known to physicists 
and engineers. Furthermore, denoting the physical dimension of an action 
by h, the Lorentz force equation (JHj shows that the absolute dimension 
of the electromagnetic field strength 2-form F is [F] = h/q = (ft, that 
is, action/charge or magnetic flux eft, see 2 . Then by ©, the absolute 
dimensions of the electric and magnetic fields E and B are [E] = <ft/t 
and [B] = (ft, respectively, and the relative dimensions [E a ] = (ft/(t£) and 
[B ab ]=(ft/f. 

Our analysis does not make use of the metric. Moreover, it is generally 
covariant and as such valid in particular in GR and SR. Furthermore it is 
valid on a spacetime manifold of arbitrary dimension. And, on top of that, 
our considerations do not depend on any particular choice of the system of 
physical units. Whatever your favorite system of units may be, our results 
will apply to it. In short: Our dimensional analysis so far is premetric, 
generally covariant, dimensionally independent, and valid for any system 
of units. 

Quite remarkably, in 1+2 spacetime dimensions we can develop the com- 
plete premetric phenomenological theory of the quantum Hall effect (QHE) 
tha t applies to the an electron gas system confined to a 2-dimensional plane 
|16j . The crucial observation is that in 1 + 2 dimensions, the current is a 
2-form, and the natural constitutive relation, assuming isotropy, is then 

J = -a u F (24) 

From the Maxwell equations © and ()1U|) we find that o"h is constant, and 
the above dimensional analysis shows that it has the dimension of conduc- 
tance. The additional microscopic study ^S] supports these conclusions 
and ultimately identifies 1/<th with the von Klitzing constant Rk- 

Since such a phenomenological scheme is premetric, the gravitational 
field neither as metric nor as linear connection ever shows up in this the- 
ory. This leads to the prediction that the QHE does not "feel" gravity. In 
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particular, this means that the linear Hall resistance should remain con- 
stant in an arbitrary noninertial frame and in any gravitational field. An 
experimental verification of this prediction is highly desirable. 

In contrast to the (1 + 2)-dimensional Hall electrodynamics, the elec- 
trodynamical theory in 1 + 3 spacetime dimensions normally incorporates 
the metric of a flat or a curved spacetime via the constitutive relation 
H = H(F) between the excitation and the field strength. In particular, the 
standard Maxwell-Lorentz electrodynamics arises when we assume 

H = A *F . (25) 

Here the (4-dimensional) Hodge star * is defined by the spacetime metric 
gij which in Cartesian coordinates reads = diag(c 2 , — 1, — 1, — 1). The 
above dimensional analysis shows that £Iq := 1/Ao has the dimension of 
resistance. Usually, it is called the vac uum impedance. 

The Maxwell eq uati ons (J3J) and (|10|) , together with the Maxwell-Lorentz 
spacetime relation (|25(l . constitute the foundations of classical electrody- 
namics: 

d*F = Q J, dF = 0. (26) 

These laws, in the classical domain, are assumed to be of universal valid- 
ity. Only if vacuum polarization effects of quantum electrodynamics are 
taken care of or if hypothetical nonlocal terms emerge due to huge acceler- 
ations, the spacetime relation H = H(F) can pick up corrections yielding 
a nonlinear law. 

We would like to stress that the presence of two fundamental constants 
in the Maxwell-Lorentz theory, namely c and Oo, is, in our opinion, very 
much underestimated. A simple application can be given in the context of 
of the possible variation of fundamental physical constants. The possibility 
of time and space variations of the fundamental constants is discussed in 
the literatu re both from an experimental and a theoretical point of view, see 
(451 1461 147] , for example. Of particular interest are certain indications that 
the fine structure constan t may s lowly change on a cosmological time scale. 
Some authors, see Peres |481 149j . for example, related some experimental 
evidence of the variability of the fine structure "constant" to the change of 
the speed of light c = c(t,x a ). However, a closer look at the definition of 
the fine structure constant 

a = ^— = = — (27) 

f 2e ch 2hA 2R K { ' 

shows, see Tobar [2H] and |29| . that it is explicitly given in terms of the ratio 
of two resistances — vacuum impedance VLq and von Klitzing constant i?x 
(the quantum Hall resistance). Note that the speed of light c disappeared 
completely! Tobar demonstrated explicitly that (27) doesn't depend on 
the system of units. This is inbuilt in our formalism right from its axiomatic 
beginnings 0. 

In other words, the formula (|27|) demonstrates that of the two funda- 
mental constants of electrodynamics, which appear naturally in Maxwell- 
Lorentz electrodynamics (see the previous paragraph), it is the vacuum 
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impedance that enters the fine structure constant and not the speed of 
light. Accordingly, a variation of the fine structure constant «f = Q!f(t) 
would force us to conclude that most probably Ao = Ao(i). A similar for- 
malism was actually developed by Bekenstein |5U] . although he inclines to 
a different physical interpretation of a variable electron charge e. Since a 
variable e and/or h would yield to charge and/or flux violation, at least one 
Maxwell equation had to be given up. Since we consider this undesirable, 
we opt for Aq = Ao(t). 



7 Axiom 5: Local and linear spacetime relation 

A local and linear spacetime relation appears to be a reasonable physical 
assumption in the general setting of the axiomatic approach. Then, the 
electromagnetic excitation and the field strength are related by the local 
and linear constitutive law 

Hij = l - Kij kl F kl . (28) 

The constitutive tensor k has 36 independent components. One can decom- 
pose this object into its irreducible pieces. Obviously, within the premetric 
framework, contraction is the only tool for such a decomposition. Following 
Post nm, we can define the contracted tensor of type [}] 

Ki k := K U kl , (29) 

with 16 independent components. The second contraction yields the pseudo- 
scalar function 

K := n k k = K H kl . (30) 

The traceless piece 

fa k := - X -Kb k (31) 

has 15 independent components. These pieces can now be subtracted out 
from the original constitutive tensor. Then, 

K ..kl _ (l) fci , (2) kl , (3) kl ( o 2 \ 
= «« ii w + 2^ fc 4] + ^%- (33) 

By construction, ' ' the totally traceless part of the constitutive 

tensor: 

( W = 0. (34) 

Thus, we split k according to 36 = 20 + 15 + 1, and the [|] tensor W Ki kl 
is subject to the 16 constraints (|34^1 and carries 20 = 36 — 16 components. 

One may call ^'Kij kl the principal or the metric-dilaton part of the 
constitutive law. Without such a term, electromagnetic waves are ruled 
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out, see jHISHl- We further identify the two other irreducible parts with 
a skewon and an axion field, respectively. Conventionally, the skewon and 
the axion fields are introduced by 

& j = -~M « = ^ K - (35) 

The standard Maxwell-Lorentz electrodynamics (|25|) arises when both 
skewon and axion vanish, whereas 

«/c^ = Ao%/<. (36) 

Here 7/ yw := yf^ghjki and rjij kl = rjij mn g mk g nl . 

Along with the original K-tensor, it is convenient to introduce an alter- 
native representation of the constitutive tensor, see Post 

X ijkl -=l^ mn K mn kl . (37) 
Substituting (|33|) into (|37|) . we find the decomposition 

^ijkl _ (l)y,ijkl _|_ (2)^ijkl _j_ (3)^ijkl ^ 

again with principal, skewon, and axion pieces. They are determined by 

(1) x i jk l = l e i jmn {l) Kmn kl } (39) 

(2) x ijkl = }_ e ijmn (2) Kmn kl = _ e ijm[k ^J}^ ^ 

(3) x ijH = 1 e ijmn (3) Kmn kl = ^e ijH K. (41) 

Using the S-identity and the K- identity derived in jSH], we can verify 
that is skew- symmetric under the exchange of the first and the second 
index pair, whereas ^'x ls symmetric: 

(2) x ijkl = _ (2) x Hij 5 (l) x ijkl = (l) x klij_ ( 42 ) 

Performing a (1 + 3)-decomposition of covariant electrodynamics 0, we 
can write H and F as column 6- vectors with the components built from the 
magnetic and electric excitation 3- vectors H a , T> a and electric and magn etic 
field strengths E a , B a , respectively. Then the linear spacetime relation ()28|) 
reads: 

H a \ ( C\ B ba \f-E b 
V a j V A ha V b a I \ B 



ha n » I I r>'< I ■ ( 43 ) 



Here the constitutive tensor is conveniently represented by the 6x6 matrix 

K ( C b a &ba \ IK f &ab T^a 



= I ) , x" = ( cT I ■ (44) 
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The constitutive 3x3 matrices A, B,C, T> are constructed from the compo- 
nents of the original constitutive tensor as 

A ba ._ x 0a 0b ^ B ba :=\e acd e befX cdef , (45) 

C\ := l -e bcd X cd0 \ V a b := l -e acdX ° hcd . (46) 
If we resolve with respect to x, we find the inverse formulas 

x 0a0b = A ba ) x abcd = e abe e cd f B fe , (47) 

,Oabc hcd-j\ a , abOc r abd pc / /i o \ 

The contributions of the principal, the skewon, and the axion parts to 
the above constitutive 3-matrices can be written explicitly as 

A ab = _ e ab_ e abcpO^ ^ 
^afe = pti + iabcPo C , (50) 

C a fe = 7 a fe - (h a ~ % Pc) + a5l (51) 
V a b = Ja+tfa" + (52) 

The set of the symmetric matrices e ab = e ba and = together with 
the traceless matrix 7% (i.e., 7 C C = 0) comprise the principal part W^-?*^ 
of the constitutive tensor. Usually, e ab is called permittivity tensor and 
fi~ b reciprocal permeability tensor ("impermeability" tensor), since they 
describe the polarization and the magnetization of a medium, respectively. 
The magnetoelectric cross-ter m 7 % is r elat ed to the Presnel-Fizeau effects. 
The skewon contributions in (|49|) and Q5U|) are responsible for the elec tric 
and mag netic Faraday effects, respectively, whereas skewon terms in Q51j) 
and (|52|) describe optical activity. The particular case of the spatially 

isotropic skewon, jS a b = | S b , j$o c = 0, jfi c ° = 0, was studied first by Nieves 
and Pal |H2] , who treated s as a third fundamental constant along with the 
vacuum impedance Oq and the speed of light c. 



8 Light propagation in vacuum 

Birefringence effects are usually studied in the geometric optics approxi- 
mation. It is equivalent to the Hadamard approach, where one studies the 
propagation of a discontinuity in the first derivative of the electromagnetic 
field. The basic notions are then the fields of the wave covector and the 
ray vector that encode the information about the propagation of a wave in 
a spacetime with a general constitutive relation. 

The crucial observation about the surface of discontinuity S (defined 
locally by a function <3? such that = const on S) is that across S the ge- 
ometric Hadamard conditions are satisfied for the components of the elec- 
tromagnetic field and their derivatives: [Fy] = 0, [diFjk] = q% fjk, [Hij] = 
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0, [diHjk] = q% hjk- Here qi := d-^ is the wave covector. Then using the 
Maxwell equations (J3J) and (|1U|) and the constitutive law (|28[). we find a 
system of algebraic equations for the jump functions: 

X l3kl q,fki = 0, -/,./;,, 0. (53) 

Solving the last equation in by means of fij = q%aj — qjdi, which yields 
the gauge invariance dj — > + we are able to reduce ()53jh eventually to 

X ijhl qjQkai = . (54) 

This algebraic system, with the mentioned inbuilt gauge invariance, has a 
nontrivial solution for aj only if the determinant of the matrix on the left 
hand side vanishes. After removing the gauge freedom in the determinant 
H|, we arrive at the 4-dimensionally invariant extended Fresnel equation 

Q ijkl (x)qiQjQkqi = 0, (55) 

with the fourth rank Tamm-Rubilar (TR) tensor density of weight +1 de- 
fined by 

G l3k \x) ■= j { e mnpq erstu X mnr(l X 3lpslk X l)qtu • (56) 

It is totally symmetric, G^ kl ( X ) = G^ kl \x)i an d thus has 35 independent 
components. 

Different irreducible parts of the constitutive tensor (|3*3*)) contribute 
dif fere ntly to the general Fresnel equation. A straightforward analysis 
[21l38| shows that the axion piece drops out completely from the TR-tensor, 
whereas the tw o remaining irreducible parts of the constitutive tensor con- 
tributes to (|Sf))) as follows: 

G ijkl ( X ) = Q t3k \ [l) X) + Wx™^ Prn Pn ] • (57) 

Birefringence (or double refraction) is the direct physical consequ ence 
of the Fresnel equation. In this case the quartic wave covector surface (|55|) 
reduces to the pair of second order light-cones. However, the influence of 
a skewon field on the Fresnel wave surface is qualitatively different: The 
characteristic sign of the skewon is the emergence of the specific holes in 
the quartic Fresnel surfaces that correspond to the directions in space along 
which the wave propagation is damped out completely jHHi- This effect is 
in complete agreement wi th our earlier conclusion on the dissipative nature 
of the skewon field 0E5|. 

9 No birefringence and the light cone 

It has been shown recently that taking the linear spacetime relation 
for granted, one ends up at a Riemannian lightcone provided one for bids 
birefringence in vacuum, see also |53| . For this the covariant equation (1551) 
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is expanded for the zeroth component qo of the 4- wave covector q = (qo, q a )'- 

M q$ + Miq% + M 2 q% + M 3 qo + M 4 = 0, (58) 

where the coefficients Mj are homogeneous functions of degree i in the 
spatial compone nts q a . Due to Ferrari (1545), the four solutions of the 
quartic equation (|58f) can be written as 



o(i)= y/a + JP+-n=-S, qI {2 )= ^-a// 3 + 4=-< 5 ' ( 59 ) 



00(1) = -V«+ JP- -j= ~S, g,J (2) = -y/a- JP- -t= -5, (60) 



where a, (3, 7, 6 depend on the coefficients Mj (25] . 

Vanishing birefringence means that there is only one future and only 
one past directing light cone. There are two possibilities, namely 

00(1) = 0(1(2) . 9o(i) = 9o(2) > Le -> /5 = 7 = , (61) 
0(1(1) = ?o(i) ' 9(1(2) = 9o( 2 ) ' Le -> a = 7 = . (62) 
Accordingly, the quartic wave surface in these cases reads 

[(qo-ql)(qo-q l o)} 2 = 0, (63) 
or, explicitly, dropping the square, 

2 1 M a If M ab M a M b \ . . 

9 3 Wi ■= qo + 2 ^fWa + g (4^- - -j^)q a q b = , (64) 

where M a , M ab are constructed in terms of the components of x^ k \ see [21 • 
Since this relation describes an unique light cone, it should be under- 
stood, up to a scalar factor, as an operational definition of a Riemannian 
metric. From the condition of t he e xistence of a unique solution (or from 
hyperbolicity, see Beig jSSl)) ea.l|64|) has to possess two real solutions for 
any given q a . As a consequence, the signature of the metric is Lorentzian. 
Accordingly, the signature of the metric is a consequence of the existence 
of a unique real solution of the Maxwell equations in a future causal cone 
for arbitrary sources with compact support. 

An alternati ve t reatment of the correlation between the linear consti- 
tutive relations (|28[) and the signature is given by two of us in [23] • In 
order to provide a physical interpretation of the 4-dimensional quantities, 
we construct their (l+3)-decompositions with a number of free sign factors. 
For the current, we take 

J = ir j A da + i s p (65) 
and for the electromagnetic field 

H = hi H A da + h s V , F = f T E A da + f s B . (66) 
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We introduced here the Time and Space factors ij, is, h T , hsfj, , fs = ±1. 
For all possible signatures of the 4-dimensional metric, we obtain the ex- 
pressions for the electric and magnetic energy densities that correspond to 
1)22(1 . The metric of a Lorentzian type turns out to be related to a positive 
electromagnetic energy density. This result does not depend on the values 
of the sign factors. 

We analyzed the (1 + 3)-decompositions of the field equations, and we 
derived which sign factors are conventional and which do depend on the 
signature. We find that the electric charge has two possible signs for all 
signatures. For all signatures, we determine the features of the interac- 
tions between charges and between currents. Only for a metric with a 
Lorentzian signature we have ordinary Maxwell-Lorentz electrodynamics. 
In particular, it yields attraction between opposite charges and repulsion 
between charges of the same sign: this is Dufay's law. Also the magnetic 
force has a correct sign. In particular, it is responsible for the pulling of 
a ferromagnetic core into a solenoid independently of the direction of the 
current, in accordance with Lenz's rule. That is, we show in the metric- 
free approach that positive electromagnetic energy density together with 
the correct signs in Dufay's and Lenz's rules correspond to a Lorentzian 
signature of the metric. 

10 Axion electrodynamics and the CFJ model 

When the skewon is trivial, jS^ = 0, the structure of the electromagnetic 
theo ry simplifies greatly. If, furthermore, the principal part has the form 
(|36j) . we end up with the linear constitutive law H = Ao*F + aF. This 
framew ork is called axio n (Maxwell-Lorentz) electrodynamics, see Ni [541 
153 OH and Wilczek j3H, e.g.: 



It is as if the current J picked up an additional piece depending on the 
gradient of the axion field. In tensor calculus, we have for the inho- 
mogeneous Maxwell equation \ Q dj(y/=g F ij ) + e ijk \d j a)F k i = J\ with 




As a degenerate special case, we can also consider the pure ("stand- 
alone") axion field with ^Ki kl = ^Kij kl = 0. Then, 




(67) 




H = aF 



or 




(68) 



and the Maxwell equations read 



(da) A F = J 



and 



dF = 0. 



(69) 



This is a special case of axion electrodynamics, namely (|67() with Ao = 0. 
Historically, the first person to discuss (and to reject) a constant pure axion 
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field was Schrodinge r I57j . p. 25, penultim ate para graph, and, as a non- 
constant field, Dic ke [58| . The framework (|68 |1 .(|69 |) in fact corresponds to 
Tellegen's gyrator 59, 60 and to Lindell & Sihvola's perfect electromagnetic 
conductor (PEMC jpHp . 

A further specialization of the axion electrodynamics is possible when 

the covector v := da = i>i dx l has constant components. In the cosmological 
context, this yields the Lorentz-violating CFJ- electrodynamics , see Carroll, 
Field, and Jackiw 

Even though in the extended Fresnel equation (|55|) the birefringence 
effect is independent of the axion field, such a behavior was discovered in 
|24j . In order to resolve this problem, we consider the standard wave ansatz 

F = fe^, (70) 

where i is the imaginary unit and <p = ip(x k ), while / is a constant 2-form. 
We denote the wave covector as q = dip = qi dx l . For the ansatz ()70j) , the 
components of the excitation 2-form for the CFJ case become 

with H kl := e klmn H mn /2. In contrast to the Hadamard method, the am- 
plitude of H kl is not a constant, even if f mn , the amplitude of F mn , is a 
constant. Substituting (|7U|) and (f?Tj) into the field equations and putting 
the current J to zero, we obtain a system of 8 linear equations 

e ijkl qjf k i = , (>'%• - i(d jX ijkl )) hi = (72) 

for 6 independent variables fki- For the special case appearing in the CFJ- 
model, the constitutive tensor involves a principal part ^x^ kl like in con- 
ventional vacuum electrodynamics and the variable axion part ( 3 )^ w = 
a(x m )e l ^ kl . Following the procedure given in [2], we rewrite this equation 
in the covariant form 

Q i3hl (x) mikdl - X i3hl (dia)(d ia ) q jqk = . (73) 

Substituting here the CFJ constitutive tensor, we obtain (f, = diet) 

{m l ?-{m l ? + {w l ){w 3 ) = Q- (74) 

Finally, if we choose a = fir, with \i as a constant and r as proper time, 
we find Vi = (/i, 0,0,0) and, together with qi = (u>,q), 

(^ 2 -q 2 ) 2 -^q 2 = 0, (75) 

which coincides with the CFJ dispersion law |24j . 

In fact we have here two different types of the birefringence effects: (i) 
The premetric birefringence is generated by the algebraic structure of the 
constitutive tensor, (ii) The CFJ birefringence is generated by derivatives 
of the constitutive tensor. 
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11 Axiom 6: Splitting of the electric current 

In order to discuss the electrodynamics of continuous media, we need some 
further input. The crucial point is as follows: The total current density 
is the sum of the two contributions originating "from the inside" of the 
medium (which is interpreted as a bound or material charge) and "from 
the outside" (which is free or external charge): 

j = jmat + je xt ( ? g) 

Here, the bound electric current inside matter is denoted by mat and the 
external current by ext. The same notational scheme is also applied to the 
excitation H, so we have H mat and H ext . 

Bound charges and bound currents are inherent characteristics of matter 
determined by the medium itself. They only emerge inside the medium. 
In contrast, free charges and free currents in general appear outside and 
inside matter. They can be prepared for a specific purpose by a suitable 
experimental arrangement (a beam of charged particles, say, and scatter 
them at the medium), or we could study the reaction of a medium in 
response to a prescribed configuration of charges and currents, J ext . 

Furthermore, we assume that the charge bound by matter fulfills the 
usual charge conservation law separately: 

dJ raat =0 ( 7? ) 

We call ((75]) together with ((7T|) Axiom 6. It specifies the properties of 
the classical material medium. In vie w o f the relation dJ = 0, resulting 
from the first axiom, the assumption (|77|) means that there is no physical 
exchange (or conversion) between the bound and free charges. Although 
the sixth axiom certainly does not exhaust all possible types of material 
media, it is valid for a sufficiently wide class of media. 

Analogously to the Maxwell equation Q, which was derived from the 
conservation law we introduce by means of (|77|) the excitation H ma,t as 
a "potential" for the bound current: 

J mat = dH mat . (78) 

The (1 + 3)-decomposition, following the pattern of Q, yields 

^mat = _ w mat A ^ + p mat_ ^ 

The conventional names for these newly introduced excitations are polar- 
ization 2-form P and magnetization 1-form M, i.e., 

pmat = _ p 5 ^mat = M _ (80) 

The minus sign is chosen in accordance w ith the usual behavior of para- 
magnetic matter. Then, by using this in (|77l) . we find, in analogy to the 
inhomogeneous Maxwell equations ©, 

-dP = p mat , dM + P = j mat . (81) 
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The identifications (|80|) are only true up to an exact form. However, the 
uniqueness is guaranteed if we require 2? mat = for E = and 7^ mat = 
for B = 0, see 0. 

In order to finalize the scheme, we define the external excitation 




H:=H-H mat ( (82) 
y$j:=H — W mat = H-M 

The external excitation M = (S),D) can be understood as an auxiliary 
quan tity. When we differentiate (|82j) and elimina te (i ff and dH m by © 
and (JToJ), respectively, we find, making use of (JTHJ), the inhomogeneous 
Maxwell equation for matter: 

dH = J cxt I d® = p™\ 
\dfi-®=f xt . 

In Maxwell-Lorentz ele ctrodynamics, we obtain from (|82|) and the universal 
spacetime relation (|23|) the expressions 

£ = £ *E + P(E,B), (84) 
£ = — *B — M(B, E) . (85) 

The polarization P(i£, P) is a functional of the electromagnetic field 
strengths E and B. In general, it can depend also on the temperature 
T and possibly on other thermodynamic variables specifying the material 
continuum under conside rati on; similar remarks apply to the magnetization 
M(B,E). The system (j83|l l ooks similar to the Maxwell equations (J5J). 
However, the equations in (J83|) refer only to the external fields and sources. 
We stress that the homogeneous Maxwell equation remains valid in its 
original form. 



12 Coupling of electrodynamics to gravity 

The coupling between the electromagnetic and the gravitational field is an 
age-old problem. It is already related to the first observable prediction 
of GR about the bending of light rays of stars in the gravitational field 
of the Sun. The electromagnetic and gravitational effects are of rather 
different orders of magnitude. However, the increasing precision of modern 
experimental techniques gives rise to the hope that the appropriate form 
of the coupling can soon be determined. 

In particular, we have in this context two independent but closely re- 
lated problems: 

(i) How does the gravitational field of a massive source depend on its 
electric charge? 



20 



F.W. Hehl, Y. Itin, and Yu.N. Obukhov 



(ii) How does the electromagnetic field of a charged massive source change 
when the gravity is "switched on"? 

In most cases, the coupling between two fields can be represented by a spe- 
cific term in the total action. This term has to respect the diffeomorphism 
invariance related to gravity as well as the gauge invariance of electrody- 
namics. Moreover, it is reasonable to require the main facts of both theories 
(the conservation of energy-momentum, of electric charge, and of magnetic 
flux) to be preserved in the modified model. Although, coordinate and 
gauge invariance strongly restrict the variety of admissible coupling terms, 
we have still an infinity set of possibilities. Indeed, we can always construct 
a polynomial of a chosen coupling term. This new term can also serve as an 
admissible additional piece of the action. Hence we will restrict ourselves 
to t he a dmissible parity conserving coupling terms of the lowest order, see 
also Ql]. 

Our main result will be that all such lowest order modifications of 
the standard Einstein-Maxwell system are completely embedded in the ax- 
iomatic approach to electrodynamics. 

12.1 Coupling of electrodynamics to Einsteinian gravity 

In the framework of GR, the coupling between the electromagnetic field 
and gravity is managed by the electromagnetic action itself 

S(g,F) = -^ J*FAF = -^J F^^-gotx. (86) 

Here g = gij dx l ®dx^ = o a p r d a ®'d^ is the metric tensor, whereas the 2- form 
F = Fijdx 1 A dxi /2 = F a p-d a A #^/2 is the electromagnetic field strength. 
The typical form of the coupling Lagrangian (|86[) is L(g,F) ~ (g 2 ■ F 2 ). 
Here and later, we use the notation ( _ • _ ) for a summation that is evaluated 
by contracting th e in dices. 

If one adds to (|86|) the actions of the gravitational and the matter fields, 
then variation with respect to the metric yields the Einstein field equation 
(without cosmological constant) 1 

Ri Cij - \ R 9ij = ^-(^ m % + (^%) , (87) 

with the Ricci tensor Ricjj := Rkij k and the curvature scalar R := g^Hicij. 
The electromagnetic and the material energy-momentum tensors act as 
sources of the gravitational field. 

1 In exterior calculus, we vary with respect to the orthonormal coframe i3 a and find 
f^^Ai^ 7 = ^( (em) E Q + (mat) E Q ). Herer/ a/37 = *(ti a A# T ), R " 3 is the curvature 
2-form of the Riemannian spacetime, and G Newton's gravitational constant. An energy- 
momentum 3-form translates into an energy-momentum tensor according to E a = T^rip, 
with r/0 = *df}. 
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The Reissner-Nordstrom solution of the Einstein-Maxwell equations 
ds 2 = / 1 _2^ + ^ dt 2_/ 1 _2^ + ^ V-rW, (88) 

with m = GM/c 2 and q 2 = GQ 2 /(47reoc 4 ), describes a mass M with an 
electric charge Q. The electromagnetic field strength F is involved in (|86[) 
only algebraically, that is, without derivatives. Hence the inhomogeneous 
Maxwell equation can only be derived by a variation taken with respect to 
the potential A. 

Let us turn then to a discussion of the potential. In the conventional 
textbook approach, see Landau and Lifshitz jSJ and Misner, Thorne, and 
Wheeler |62| . the components Ai of the electromagnetic potential 1-form 
A = Ai dx l are taken as field variables. The recipe for the "minimal" 
tran sition from SR to GR is then the "comma goes to semicolon rule" 
|62| : the comma "," denotes the partial, the semicolon ";" the covariant 
derivative. For the components Fu of the field strength F we find then 

F%j — — > (^9) 

since the Levi-Civita connection is symmetric. The relation between the 
components Fij and A{ is then eventually recognized as independent of the 
metric structure. Needless to say that such a procedure is highly coordinate 
dependent. 

In contrast, in our axiomatic approach the potential 1-form A, as (coor- 
dinate independent) geometrical object, is the field variable of our choice. 
From magnetic flux conservation we find directly 

F = dA , (90) 

without messing around with covariant derivatives. In this coordinate- 
independent way we recognize right away that there is no chance for the 
metric to intervene in F = dA. Clearly, a coordinate independent proce- 
dure is to be preferred against a coordinate dependent one. 

The conventional story then continues as follows: We "semicolonize" 
the Maxwell equations in order to make them fit for survival in the curved 
pseudo-Riemannian spacetime of GR, 

F^.j = n J i - 

In the first equation the metric tensor drops out, again because of the 
symmetry of the Levi-Civita connection: 

Fij,k + Fjk,i + Fki,j = . (92) 

It is not sensitive to the gravity field of the source. The second equation 
can be rewritten as 



(93) 
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Hence the gravitational field is involved in the inhomogeneous Maxwell 
equation only via the metric dependent expression within the parenthesis 
and via the determinant of the metric tensor. Incidentally, since for the 
Reissner-Nordstrom solution the determinant is the same as for a flat mani- 
fold, the electromagnetic field is not sensitive to the mass of the point 
charge. 

In our axiomatic approach, see (|26j). we arrive directly at 

f dF = , 
\d*F = n J, 

without discussing covariant derivatives. The Hodge star is known to de- 
pend only on the conformal part of the metric, see Frankel jlj. If we 
identify the currents according to J % = \/—ge l ^ kl Jjki/6, here J = Jijkdx 1 A 

dxi A dx k /6 is the cur r ent 3-form , th en the endresults of both procedures 
coincide; compare pZ |) .l|9"5 |) with IjSijl. 

The Maxwell equations in tensor notation, if altered in the presence of 
the gravitational fiel d, c an be rearranged by a change of the basic notation. 
In accordance with (|25|). the electromagnetic excitation is linearly related 
to the field strength: 

H ij = ~ x iiU F k i . (95) 
The constitutive tensor reads 

x m (g) = A V^(g lh g> 1 - g ll g 3k ) • (96) 

Then the electromagnetic action (|86p takes the form 

S(g, F) = — i J EijfrWx = -~jFAH, (97) 

whereas the field equations in tensor analytical and in exterior form nota- 
tions become, respectively, 

J F ijJk + F jk ,i + F fci j = , or | dF = , 

|^i ii = J% " \dH = J. 

Consequently, in the Einstein-Maxwell syst em, the coupling is completely 
determined by the constitutive tensor (|96|) . This coupling is referred to 
as minimal coupling. Technically, the minimal coupling is provided in 
the conventional coordinate dependent textbook approach by the substi- 
tution of the partial by the covariant derivatives (taken with respect to the 
Levi-Civita connection) and in the axiomatic exterior calculus approach by 
Hodge staring the field strength F. 

Summing up we can state that t he m ini mal coupli ng b etween electrody- 
namics and gravity is provided by (|2^)l [or (|9"5)) with (|95|)] . The knowledge 
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about the Riemannian structure of spacetime is fed into the Hodge star 
operator via the metric. In SR a flat metric appears, in GR a curved one. 

Nonminimal coupling 

If for some reason we believe that the minimal coupling between elec- 
trodynamics and gravity i s in sufficient and has to be generalized, then we 
have to modify the action (|86[) by adding so-called nonminimal terms. Still, 
we have to respect diffeomorphism and gauge invariance. Riemannian ge- 
ometry is characterized by the curvature tensor Rijk and its traces Rib- 
and R. Terms proportional to the electromagnetic potential are forbidden 
because of gauge invariance. The term linear in F vanishes as a product of 
a symmetric and an antisymmetric tensors. Consequently the lowest order 
nonminimal Lagrangian reads 

L(F, R) = a 1 F ij F kl R i ^ kl + a 2 F ik F^m C i j + a 3 i^F ij R . (99) 

Let us recall several arguments in favor of a such non-minimal couplings: 

(i) For the values a\ = 03 = — 4ct2, eq.(l9"9l can be recovered from the 
5-dimensional Gauss-Bonnet action 155] . 

(ii) Another special set of parameters in Q99JI 

e 2 A 2 13e 2 Ag e 2 A c 2 

ai = "7^h' a2= 7^h' a3 = -288^h' (1 ° 0) 

were de rive d in the one-loop approximation of quantum electrodynamics 
(QED) |67U68j . Here e is the elementary charge, h = 2tvH Planck's constant, 
and A c = h/(mc) the Compton wavelength of the electron. 

(iii) The Lagrangi an ()99|) was used as a basis for models with a variable 
speed of light [59l 1701 T7T] . Such an effect is present in all models, except 
for the cases a± = a 2 = 0. 

The constitutive tensor corresponding to the Lagrangian L(F, g)+L(F, R) 
depends on the metric and the curvature tensor and its contractions: 

x ijkl (g, R) = (1 - 4a 3 WV - gV) - 4<*iR ijkl 

-a 2 (g ik Ric jl - g il Ric jk + g jl Ric ik - g jk Ric il ) . (101) 

Since (jlUljl is derived from a Lagrangian, i ts s kewon part (|40[) is zero. Also 
the completely antisymmetric axion part (|41|) vanishes due to the symme- 
tries of Rijki- Thus, in (jlOlj) only the principal part ^^ kl {g^R) is left 
over. Substituting (|1U1|) into the 3x3 constitutive matrices (|45 |) .(|46j ) . we 
obtain A = A T , B = B T , C = T> T . The a 1 and a 2 terms yield birefringence 
in general. 



24 



F.W. Hehl, Y. Itin, and Yu.N. Obukhov 



12.2 Maxwell's field coupled to Einstein-Cartan gravity 

In Einstein-Cartan gravity, see Blagojevic |72| . e.g., the spacetime geometry 
is of the Riemann-Cartan type and as such on two fundamental structures 
- the metric tensor gij and the connection Tij k . These quantities satisfy 
the metricity condition 

9ij;k = 9ij,k — ^kidlj — Fk/gil = . (102) 

If the connection is not symmetric, then the torsion tensor Ti k is in- 
volved. In holonomic coordinates, it is defined by T{ k = 2T^jj k . In exterior 

calculus, we have the 2-form T a = ZM a = Tij a dx % A dx 3 /2, with the exte- 
rior covariant derivative D. Accordingly, a possibl e interaction of torsion 
with the electromagnetic field is of special interest [7J3 EU ESI 123 IUJ EH 

The minimal way of coupling the electromagnetic field with Einstein- 
Cartan gravity works the same way as with Einstein gravity. Since neither 
the exterior derivative d nor the H odg e star * can feel torsion directly, we 
have again the minimal equations (|26|) : 

d*F = n J, dF = 0. (103) 

They will do the job. 

The conventional coordinate dependent textbook approach led to nu- 
merous misunderstandings. If one applies the semicolon rule to the special 
relativistic formula Fij = Aji — Aij, then one arrives at 

-Aj-i Ai-j = Ajj Aij T{j Ah = Fij Tij A^ . (104) 

This expression is diffeomorphism invariant. However, it is not gauge in- 
variant under Ai — > Ai + di<p. Consequently, it has to be rejected. As we 
saw in the last paragraph, the axiomatic exterior calculus approach leaves 
the definition F = dA intact, and there is no need, nor is it allowed, to 
entertain in covariant derivatives in this context. 

Accordingly, the free Maxwell Lagrangian is added to the gravitational 
Lagrangian with the unchanged field strength. Nevertheless, in an exact 
solution of a gravity-Maxwell system, the torsion may depend on the elec- 
trical charge, as is exemplified by the Reissner-Nor dstr om solution with 
torsion in the framework of a Poincare gauge theory |79| . In other words, 
torsion is influenced by the electric charge indirectly via the field equa- 
tions. What in the Lagrangian looks like "no interaction at all," still yields 
an effective "minimal" interaction. 

Nonminimal coupling 

We are now looking for a complete family of nonminimally coupled 
Maxwell-torsion Lagrangians that are diffeomorphism and gauge invariant. 
Because of gauge invariance, the potential Ai must not appear in the action. 
Moreover, all expressions have to contain an even number of indices in order, 
if contracted, to yield a scalar. The expressions linear in the field strength 
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are of the form (T • (g 2 • T 2 )). Although such terms yield a family of 
diffeomorphism and gauge invariant Lagrangians, they have to be rejected 
from a physical point of view. Indeed, on the level of the field equation, 
the Lagrangians linear in F admit the existence of a global electromagnetic 
field even without charges. Such a modification of classical electrodynamics 
seems to be unwarranted. 

Hence the lowest order addenda to the Lagrangian are quadratic in F; L j 

and quadratic in torsion Tij k , i.e., of the typical form F 2 ■ g 3 • T 2 . The 
contractions can be performed in the following three ways: 

(i) All indices of the F-pair and of the T-pair are contracted separately, 



(F 2 .g 2 )(g.T 2 ), 



or, explicitly: 

WL(F,T) 
P>L(F,T) 



JT. . Tpij T , T mnk 
1 ij A - 1 mnk- L j 

JT. ■ JTV T nrpmk 
it 1 - J -mn ± k 



(2) 



L(F, T) — FijF lj T mnk T 



nkm 



(105) 



(ii) Two free indices of the .F-pair are contracted with two free indices of 
the T-pair, i.e., 



(f ■ g) ■ {g ■ t ) . 



JP. jr. m Tinjrp k 
1 ira 1 'J ± ± nk 



Such terms are 
WL(F,T) = 
V>L(F,T) = 
WL(F,T) = 

(iii) The F-pair and the T-pair have four free indices: 



JP. jp.m rpiknrpj 
1 ira 1 j ± -L nk ■ 

m rrtik 



jr. jr. m T™,Ti r 

1 im 1 j 1 k ± 



^L(F,T)=F m F j m T ikn Thn. 

m rpiknrp 
1 J-k 

(106) 



^L(F,T) = F im F J m r kn T k J 



(F 2 ).(g 3 -T 2 ), 



WL(F,T) = 

(11) L(F,T) = 

WL(F,T) = 

(15) T(F,T) = 

WL(F,T) = 



kl 



jr. . jr. , T im T 3 

1 ij 1 kl m± 

F tJ F kl T m ^T mkl , 

1 ij 1 kl ± m ± i 



FijFkl T 

FijFkl T 



mik 



ikm 



-L 1 



WL(F,T) 
(12) L(F,T) : 
^L(F,T) . 
(16) T(F,T) 



Summing up, the general torsion Lagrangian reads 

17 



- JT JT n^ini rpklj 
' ij kl 1 m-L j 

F l3 F kl T m lk T ml \ 
F ■ Fui T mij T kl 

1 I] 1 kl ± 1 mi 

--F,F kl T^ m T kl m , 

(107) 



(0 



L(F,T), 



(108) 



=i 
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with the dimensionless constants Bj. 

A special Lagrangian of type (jlU8j) . namely 

*(T a A F) T a A F = i ( - 4 (5) L + (16) l) *1 (109) 

was considered recently [171 118j in the context of a cosmological test of 
Einstein's equivalence principle. 

Since we start with a modified Lagrangian, a skewon piece does not 
occur in the corresponding constitutive tensor. Thus 

x m^ T )=i 2 [Y^ (T mnp T qrs ) } MW = x ijk \g, T) + x ijk \g, T) , (110) 

where x* J (<7, T) is the modification of the principal part, while the axion 
part takes the form 

X ij %,T) =i 2 [Y, (T m n P T qrs )] [tjkl] . (HI) 

Birefringence induced by torsion. For the Lagrangians (|105[) - (|107|) 

we find the following effects on x^ kl and on the light cone (we put I = 1): 

(i) The Lagrangians (|105j) yield 

r ikl (9, T) = S(g lh gi l - g il g^) = 2^V ]i , (112) 

where S is a scalar function quadratic in torsion. Therefore, the Tamm- 
Rubilar tensor density changes only by a conformal factor and the light 
cone is preserved. 

(ii) For (EH) , we introduce the abbreviation S ij = S ji := [g 2 • T 2 ] . The 
axion field is absent. Thus, 

X ijkl (g,T) = g^ k S lM . (113) 
The 3x3 constitutive matrices A,B,C,V read: 

A ab (g,T) := x° b ° a (g, T) = \(g 00 S ab + g ab S 00 ) , 
B ab (g, T) := - A e bcd e aef x° def (g, T) = \ {g ab S c c - S ab ) , 
C\(g,T) := ^e bcd x cd0a (g,T) = \e a bc S 0c . (114) 
These matrices, with C a b (g,T) = D b a (g,T), obey 

A(g,T) =A T (g,T), B(g,T) = B T (g,T), C(g,T) = D T (g,T) (115) 
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(T = transposed). These relations, together with the closure condition for 
^ (g,T), guarantee the uniqueness of the light cone 0. Thus birefrin- 
gence does not emerges in this group of models either. 

(hi) For the Lagrangians from the t hird group l|107j) . birefringence is a 
generic proper ty. This was shown for (|1U9|) in the case of spherically sym- 
metric torsion [2U] . 

Axion field induced by torsion. The axion, that is, a pseudo-scalar 
field, was extensively studied in different contexts of field theory. Its quan- 
tized version is believed to provide a solution to the strong CP problem 
of quantum chromodynamics (QCD). The emergence of axions is a gen- 
eral phenomenon in superstring theory. The axion as a classical field also 
appears in various discussions of the equivalence principle in gravitational 
physics and in inflationary mod els. In the con text discussed in this paper, 
the axion field is introduced in (|41|) with 1)35(1 9 as one irreducible pi ece o f 
the electromagnetic constitutive tensor of spacetime. According to 
we may imagine that in the present case the axion field is induced by the 
nonminimal coupling of the Maxwell field to the torsion of spacetime. Ex- 
plicitly, we find the following: 

(i) The Lagrangians (|l(J5j) yield 

X ijkl (g,T)=t 2 gl ik gM(T.T) = Q. (116) 

(ii) The Lagrangians (|106j) yield 

X ijkl (g, T) = i 2 g^\T ■ Tf® = . (117) 

(hi) For those of the third group (|1U7|) . the general form of the axion field 
a = e ijklX ijkl (9,T)/4l is 

ol = fe ijkl (a{F^ m T kl m + a 2 T^ k T lm m ) , (118) 

where ai,a 2 are free dimensionless par amet ers, which are linear combina- 
tions of the /3's. For the special case (|109|) . the axion is extracted from 
(|118|) by putting a\ = 1, a 2 = 0. 

If in vacuum an axion field a emerges, then the Lagrangian picks up 
an additional piece ~ a F A F, see 0. Accordingly, the Maxwell equations 
are those displayed in (|67|) . Only a non-constant axion field contributes. 
The coupling of the Maxwell field to a non-constant axion field, in the case 
of a plane electromagnetic wave, amounts to a rotation of the polarization 
vector of the wave, see [HO], i.e., the axion field induces an optical activity. 



12.3 Maxwell's field coupled to metric-affine gravity 

Metric-affine gravity (MAG) is based on a s pac etime geometry with com- 
pletely independent metric and connection jHJ. Thus, in addition to the 
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torsion tensor, the nonmetricity tensor emerges 



Q 



9ij,k 



-9ij,k + ^kidlj + ^kj Oil ■ 



(119) 



In exterior calculus, we have the nonmetricity 1-form Q a g := —Dg a p, with 
the decomposition Q a p = Qi a pdx l . In this framework, we can look for 
a possible interaction of the nonmetricity with the electromagnetic field. 
As in Einstein and Einstein-Cartan gravity, the minimal coupling remains 
untouched, that is, we have again \od*F = J, dF = 0. 

In nonminimal coupling, however, there emerge two additional types of 
terms, namely F 2 ■ g 5 ■ Q 2 and F 2 ■ g 4 ■ Q - T. These coupling terms can be 
expanded in the following ways: 

(i) All indices of the F-pair and of the Q-pair (or of the TQ-pair) are 
contracted separately: 



(F 2 ■ g 2 )(g 3 ■ Q 2 



or 



(F 2 .g 2 )(g 2 .Q.T). 



Examples of such terms can be easily constructed: 



WL(F,Q)=F ij F i iQ mnk Q 



rank 



(1) 



L(F, Q, T) — FyF^ Q mn kT 



mnk 



(120) 



The corresponding constitutive tensors are of the form (S is a scalar) 

x ljkl (g,Q) = (g tk g jl -g u g' k )S. (m) 

Certainly the axion field and the birefringence effect are absent in this group 
of models. 

(ii) Two free indices of the F-pair are contracted with two free indices of 
the Q-pair (or TQ-pair), i.e., 



(F 2 ■ g) ■ ■ Q 2 



or 



{F ■ g) ■ (g -Q -T). 



Examples of such Lagrangians are 



(2) 



F(F, Q) — FikFj k Q m m nQ 



(2), 



L(F, T, Q) = F ik Fj T mn n Q tm: > . 

(122) 

The corresponding constitutive tensor is of the form 



(123) 



where = (g ■ Q ) lJ or S^- 7 = (g ■ Q ■ T) 13 . This tensor is necessary 
symmetric since it multiplies the symmetric tensor Fi k Fj k . Similarly to 
the non- minimal coupling to torsion (|113I115|) . the birefringence effect is 
absent in this group of models. Comparing to ()117|) we see that also the 
axion part of this constitutive tensor is zero. 
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(iii) The F-pair and the Q-pair (or TQ-pair) have four free indices: 



(F 2 )-(g 5 -Q 2 ), 



or 



(F 2 ).(g*.T.Q) 



For instance, 
^L(F,Q) = F ij F kl Q* k m Qi lm , 



(3) 



L(F,T,Q) 



F- Fi, T im O kl ■ 

1 ij ± kl - 1 m^c i • 



(124) 



Possessing an axion field and birefringence are generic properties in this 
group of models. 

In MAG, all the non- minimal coupling terms described above are of the 
same quadratic order in the connection. Therefore, additional Lagrangians 
depending on curvature, torsion, and nonmetricity have in general to be 
considered together. 



13 Discussion 

Although Maxwell's electrodynamics is a firmly established classical field 
theory, it is still open to new developments along many lines. Furthermore, 
the study of the fundamental structures of Maxwell's theory can serve as 
a natural starting point for non-abelian modifications with applications in 
high-energy physics and gravity. In particular, modern string and brane 
theories rely heavily on classical field-theoretic structures. In addition, 
there are important classical problems still awaiting for their solutions, 
such as the quantization of the electric charge, the existence of magnetic 
monopoles, and the value of the coupling constant, as well as a number of 
new open issues: (i) The physical place and the role of the premetric part- 
ners of the standard electromagnetic field, such as the axion, the skewon, 
and the dilaton fields, (ii) The appropriate form of the coupling between 
the electromagnetic and the gravitational field (also taking into account 
the possible gauge-theoretic extensions of standard GR, such as Poincare 
gravity and, more generally, metric-affine gravity), (iii) The description of 
a high-energy electromagnetic field in a material medium. 

In this lecture we gave an overview of the recent developments of the 
axiomatic premetric approach based on the conservation of electric charge 
and of magnetic flux as well as on some additional inputs such as the struc- 
ture of the energy-momentum. These facts are well established theoretically 
and tested experimentally. The field variables and the corresponding field 
equations are then straightforwardly derivable from the basic axioms when 
certain natural restrictions on the topology of spacetime are assumed. The 
resulting construction is actually a topological one, in the sense that it does 
not depend on a specific geometry of the underlying spacetime manifold. 
Only the additional notion of a constitutive relation brings in the infor- 
mation on the specific structure of the geometry of spacetime. The open 
problems, indicated above, are directly linked to the investigation of the 
properties of the constitutive relation. Going beyond the linear and local 
constitutive law is, in fact, a step which takes into account the quantum 
nature of the actual matter sources. The corresponding dynamics of the 
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electromagnetic field becomes highly nontrivial. A good recent analysis of 
the earlier nonlinear electrodynamical models can be found in I36j : specific 
results on the birefringence in such theories were obtained in |37| . 
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